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A  STATE-OF-THE-ART  REVIEW  OF  IRREGULAR  WAVE  GENERATION  AND  ANALYSIS 

BY 

Professor  Bruce  Johnson,  U.S.  Naval  Academy 
In  collaboration  with  the  Seakeeping  Committee  of  the  16th  ITTC 


INTRODUCTION 

In  1979,  a  questionnaire  on  wave 
generation  and  analysis  systems  prepared 
by  the  Hydraulics  Laboratory  of  the 
National  Research  Council  of  Canada  was 
sent  to  member  organizations  of  the  ITTC 
as  well  as  to  hydraulics  laboratories 
around  the  world.  Replies  were  received 
from  98  of  184  Hydraulics  Laboratories 
and  from  44  of  74  Institutions  which 
have  ship  model  towing  basins.  The 
results  of  this  questionnaire  were 
distributed  to  those  Institutions  which 
responded  and  were  summarized  In  a  paper 
presented  to  the  Coastal  Engineering 
Conference  In  Sydney,  Australia  In 
March  1980  (1).  Although  a  higher 
percentage  of  towing  tanks  reported 
Irregular  wave  capabilities  (67%  to 
41%),  they  also  reported  a  very  low 
utilization  of  those  capabilities  for 
ship  model  testing.  This  Is  most  likely 
related  to  the  nature  of  ship  model 
testing  as  compared  with  the  types  of 
tests  carried  out  In  hydraulics  labora¬ 
tories  and  ocean  engineering  basins.  It 
may  also  be  related  to  discussions 
concerning  the  validity  and  applica¬ 
bility  of  deterministic  vs.  stocastlc 
approaches  to  wave  generation  and 
analysis.  Consequently,  this  review  Is 
organized  according  to  the  methods  of 
Irregular  wave  generation  and  analysis. 


I.  NON -DETERMINI STIC  GENERATION  - 

STOCASTIC  ANALYSIS 

Non-determlnlstlc  wave  generation 
1$  defined  as  the  use  of  a  non-repeating 
wave  train  derived  by  some  method  from 
a  "true"  random  noise  source  such  as 
thermal  or  cosmic  noise.  Stocastlc 
analysis  refers  to  the  fact  that  proba¬ 
bilistic  methods  must  be  used  to  obtain 
the  customary  time  series  analysis 
functions  such  as  spectral  density 
functions,  probability  density  functions, 
correlation  functions,  etc.  This  Is 
because  individual  estimates  of  these 
functions  for  random  phenomena  show 
considerable  variance  with  time. 

Averaging  must  be  used  to  consider  the 
analysis  as  approaching  that  of  a 
stationary  random  process.  A  stationary 
random  process  Is  one  for  which  the 
values  of  the  time  series  analysis 
functions  are  Independent  of  a  shift 
In  the  time  origin  of  the  analysis  (1). 


Until  recently,  the  most  common 
methods  of  Irregular  wave  generation 
Involved  the  use  of  non-repeating  wave 
trains  derived  from  the  use  of  a  “true" 
random  white  noise  source  (3).  The 
output  from  the  noise  generator  can  be 
shaped  by  analogue  filters  to  compen¬ 
sate  for  the  transfer  function  of  the 
wavemaker  and  to  shape  the  desired  wave 
spectrum  (4).  The  shaping  filters  can 
be  set  using  trial  and  error  procedures 
or  can  be  programmed  by  computer  to 
achieve  the  desired  spectrum  (4).  The 
resulting  wavemaker  drive  signal  Is 
generally  stored  on  analogue  or  digital 
magnetic  tape  or  on  punched  paper  tape 
for  use  when  desired.  A  digital  computer 
can  also  be  used  to  generate  the  non¬ 
repeating  random  numbers  which  could  be 
converted  Into  white  noise,  but  It  Is 
more  common  to  simulate  true  random 
waves  by  generating  very  long  pseudo¬ 
random  (repeating)  sequences  In  the  time 
domain  using  binary  shift  register 
techniques  (5,6).  A  true  continuous 
spectrum  Is  non-repeating  by  definition 
since  the  frequency  Increment  approaches 
zero.  Pseudo-random  techniques  Involv¬ 
ing  very  long  repetition  periods  can 
approach  a  nearly  continuous  spectrum 
but  they  still  represent  a  discrete 
spectrum.  The  most  frequently  used 
pseudo-random  time  series  are  station¬ 
ary*  random  processes  when  their 
statistics  are  taken  over  the  entire 
repetition  period.  This  is  because  the 
component  spectral  density  values 
average  out  to  the  constant  coefficients 
in  a  Fourier  series  expansion.  Wind 
generated  waves  In  the  actual  ocean  are 
quite  unlikely  to  have  constant  spectral 
coefficients,  however,  since  energy  Is 
being  added  by  the  time-varying  wind 
fields,  extracted  by  wavebreaking,  and 
redistributed  by  non-linear  wave  Inter¬ 
actions.  Indeed,  the  resulting  non¬ 
stationary,  time-varying  spectral 
density  components  present  during  a 
severe  storm  greatly  increase  the  possi¬ 
bility  of  generating  extreme  waves  as 
compared  with  artificial  wave  generation 
methods  Involving  In-phase  addition  of 
constant  amplitude  components  In  pseudo¬ 
random  waves.  Thus,  there  Is  some 
question  as  to  whether  this  type  of 
pseudo-random  sequence  or  even  “true 
random  Gaussian  noise"  should  be 
considered  adequate  to  simulate  non- 
determlnlstlc  waves  used  In  testing 
stationary  models  subject  to  extreme 
wave  conditions. 


Techniques  utilizing  actual  sea 
records  have  also  been  developed  (7). 

In  this  case,  the  component  spectral 
density  values  generally  vary  with  time 
and  one  has  the  problem  of  finding  a 
record  length  long  enough  to  approach  a 
stationary  random  process  If  standard 
Gaussian  statistical  properties  are  to 
be  assumed  (2).  If  "wave  grouping" 
effects  are  to  be  investigated,  there  is 
also  the  problem  of  determining  whether 
or  not  a  given  prototype  wave  record 
contains  the  possibility  of  generating 
desired  extreme  waves  at  the  model  since 
the  location  of  the  wave  height  measur¬ 
ing  device  is  frequently  not  at  the  site 
where  the  prototype  structure  Is  to  be 
located.  Therefore,  the  wave  time 
histories  which  impact  on  the  model  may 
not  represent  what  can  be  expected  on 
the  actual  site,  even  though  they  may 
have  the  same  spectrum.  This  problem 
is  discussed  further  in  section  IV. 

Since  long  term  averaging  of 
stocastic  wave  records  is  required  to 
satisfy  Chi-squared  statistics  for 
reasonable  confidence  intervals,  this 
type  of  testing  Is  more  commonly  used 
for  stationary  models.  For  tests 
involving  moving  models,  Yamanouchi  (8) 
advises  that  each  run  should  be  analyzed 
separately  and  the  results  averaged 
rather  than  attempting  to  splice  together 
the  records  from  short  runs  to  make  a 
long  record.  The  latter  Introduces 
errors  In  the  estimation  of  the  spectrum 
as  well  as  to  the  estimation  of  the 
response  amplitude  operator. 

"Windowing"  should  also  be  used  for 
the  analysis  of  non-determi ni sti c  time 
series  data  since  discrete  Fourier 
analysis  presumes  that  the  sampled  data 
repeats  each  frame.  A  "Hanning"  or 
other  type  of  window  reduces  the  series 
to  zero  amplitude  at  the  beginning  and 
end  of  the  analysis  frame.  A  Fourier 
Integral  may  be  used  to  describe  the 
non-repeating  series  of  data  points 
which  has  no  discontinuities  at  each 
end  of  the  finite  segment  of  the  time 
history  being  analyzed  (2). 

Since  any  window  (other  than 
rectangular)  effectively  reduces  the 
energy  contained  in  an  analysis  frame, 
a  correction  to  the  appropriate  time 
series  function  must  be  applied.  There 
is  some  question  on  just  how  to  apply 
this  correction  -  whether  to  normalize 
the  window  time  series  to  the  original 
energy  level  or  to  apply  a  correction 
to  each  value  of  the  spectral  density 
estimates  or  their  averages.  The 
author  was  unable  to  find  a  satisfac¬ 
tory  discussion  of  this  problem  in  the 
typical  books  and  literature  on  time 
series  analysis. 


II.  DETERMINISTIC  GENERATION  - 

STOCASTIC  ANALYSIS 

Deterministic  wave  generation  is 
taken  to  imply  that  a  repeating  wave 
train  is  generated  which  contains 
discrete,  harmonically  related  compo¬ 
nents  whose  spectral  densities  are  fixed 
in  amplitude  but  randomized  In  phase. 
This  Is  based  on  the  fact  that  any 
repetitive  time  series  may  be  expressed 
as  a  Fourier  Series  expansion  with 
constant  coefficients.  Stocastic 
analysis  (Averaging)  is  required, 
however,  anytime  the  testing  time  is 
less  than  the  repetition  period  of  the 
wave  generation  record.*  Since  the 
repetition  of  the  wave  is  the  inverse  of 
the  spectral  spacing  of  the  discrete 
generating  frequencies,  the  longer  the 
repetition  period,  the  closer  a 
continuous  spectrum  Is  approximated. 

Time  domain  generation  of  this  type  of 
wave  includes  the  mechanical  gear 
systems  of  the  type  offered  by  Kempf 
and  Remmers,  and  electronic  sinusoidal 
summing  devices,  both  of  which  have 
predetermined  numbers  of  discrete 
components  (some  of  which  may  have  2ero 
amplitude).  The  previously  mentioned 
pseudo-random  waveform  generators 
pioneered  by  the  Wallingford  Research 
Station  and  available  from  EEL  Labora¬ 
tories  of  British  Hovercraft  Corporation, 
are  generally  hardwired  to  enable  the 
user  to  choose  the  repetition  period. 
They  may  contain  manually  adjusted  or 
programmable  shaping  filters  to 
compensate  for  the  wavemaker  transfer 
function  in  obtaining  a  desired  spectrum 
in  the  tank  (4,  5,  6). 

Frequency  domain  generation  tech¬ 
niques  based  on  the  inverse  Fourier 
Transform  are  gaining  increased  popu¬ 
larity  with  those  who  have  access  to 
an  on-line  computer,  since  the  desired 
spectrum  can  be  achieved  through 
iteration  as  demonstrated  by  Anderson 
and  Johnson  (9).  Funke  and  Mansard  (10) 
have  added  an  additional  dimension  to 
this  technique  by  combining  amplitude 
and  phase  modulation  to  produce  a 
desired  grouping  of  the  waves.  This 
Is  discussed  further  under  "extreme 
wave  generation". 


♦Footnote:  Stocastic  analysis  can  also 
be  used  to  advantage  In  regular  wave 
testing  to  compensate  for  the  deficien¬ 
cies  of  certain  types  of  wavemakers 
which  produce  beat  phenomena,  etc.  In 
this  case,  the  waves  are  treated  as  a 
narrow  band  random  process. 


III.  DETERMINISTIC  GENERATION  AND 

ANALYSIS 

The  previously  mentioned  pseudo¬ 
random  wave  generation  techniques  can 
be  deterministically  analyzed  provided 
that  the  analysis  time  1$  precisely 
equal  to  one  or  more  repetition  periods 
of  the  generation  sequence.  This  "zero 
variance"  technique  (9)  has  the  advan¬ 
tage  of  being  "ergodlc"  In  the  sense 
that  the  accuracy  of  the  analysis  is 
Independent  of  which  portion  of  the  wave 
record  1$  analyzed,  and  It  does  not  have 
to  be  averaged  or  windowed  unless  the 
model  exhibits  non-ergodlc  response. 

Since  there  Is  an  Integer  number  of 
cycles  of  each  discrete  frequency  com¬ 
ponent  In  any  analysis  frame  equal  In 
length  to  the  repetition  period, 
decisions  concerning  averaging  can  be 
based  on  the  desired  number  of  wave 
encounters  to  Insure  the  model  response 
has  stabilized.  This  Is  traded  off 
against  the  possible  Influence  of  beach 
reflections  on  what  was  supposed  to  be 
a  stationary,  ergodlc  spectrum  at  the 
test  location  In  the  tank.  If  the 
testing  time  Is  not  limited  by  a  moving 
model  or  beach  reflections.  It  would 
probably  be  better  to  double  the  number 
of  spectral  lines  by  doubling  the 
repetition  period  rather  than  by  averag¬ 
ing  two  "Identical"  records.  There  Is 
still  an  unresolved  Issue  concerning 
the  minimum  number  of  spectral  lines, 
the  maximum  allowable  spectral  line 
spacing  and  maximum  allowable  wave 
height  Increments  on  histograms  for 
various  kinds  of  model  tests  which 
utilize  discrete  spectrum  wave  genera¬ 
tion.  In  mooring  studies,  for  example, 
misleading  results  can  be  obtained  from 
both  regular  wave  and  irregular  wave 
studies  with  widely  spaced  harmonic 
components  (3).  This  Is  because  non¬ 
linear  and  asymmetric  elastic  charac¬ 
teristics  of  mooring  lines  and  fenders 
along  with  slowly  varying  wave  drift 
forces  can  produce  sub  and  super¬ 
harmonic  components  not  Identified  by 
discretely  spaced  regular  or  Irregular 
wave  frequencies.  However,  the  "zero 
variance"  or  periodic  Irregular  wave 
tests  can  supply  the  same  Information 
as  an  equivalent  number  of  regular  wave 
tests  for  each  of  the  discrete  spectral 
components. 

In  order  to  make  the  analysis 
deterministic.  It  Is  Important  to  match 
the  clocks  controlling  the  wave  genera¬ 
tion  with  clocks  associated  with  the 
data  acquisition.  At  the  Naval  Academy, 
for  example,  the  data  is  digitized  by 
an  analogue  to  digital  converter 
sampling  at  frequencies  which  are 
multiples  of  powors  of  two.  This  Is  to 
Insure  that  the  Fast  Fourier  Transform 
analysis  (whose  frame  size  normally 
must  be  a  power  of  two  unless  a  different 
radix  Is  used  In  the  algorithm)  will 


yield  "nice"  frequency  spaclngs  llke.l, 
.05,  or  .01  hertz.  The  clocking  of  the 
digitized  wavemaker  drive  signal  Is  at 
25.6  samples  per  second  (smoothed  by  an 
8  pole  analog  filter)  so  that  a 
repeating  wave  of  20  seconds  duration 
(.05  Hz  spacing)  has  a  drive  signal  512 
words  long. 


This  type  of  testing  Is  still  more 
appropriate  for  stationary  models, 
however,  where  testing  time  is  not  a 
function  of  the  model  speed  and  length 
of  the  tank.  A  more  serious  problem 
faced  by  moving  models  Is  that  the 
encounter  frequency  shift  causes  a  deter¬ 
ministically  generated  wave  appear  as  a 
stocastlc  wave  to  the  moving  model. 

This  Is  because  the  equally  spaced 
generation  frequencies  are  now  unequally 
spaced  In  the  encounter  frequency 
domain;  1 .e. , 


and 


Afe  «  (1  ♦  Af  where 

f  *  wave  frequency  In  fixed 
coordinates 

f  *  encounter  frequency 
V  *  model  velocity 
g  *  gravitational  constant 


This  latter  problem  Is  currently 
being  studied  by  the  col laborators  on 
this  report,  who  are  both  working  on 
"single  pass"  seakeeping  tests  for 
moving  ship  models  In  head  or  following 
seas.  The  goal  of  such  studies  Is  to 
reduce  the  testing  time  required  for 
seakeeping  tests  which  utilize  regular 
waves.  By  replacing  the  tests  In 
regular  waves  with  tests  In  determinis¬ 
tic  Irregular  or  transient  waves,  the 
cost  of  such  tests  can  be  reduced,  or 
the  range  of  parameters  Investigated 
(such  as  non-linearity)  can  be  Increased 
for  the  same  amount  of  tank  time. 
Admittedly,  such  tests  cannot  handle 
deck  wetness  and  slamming  probabilities 
directly,  but  the  resulting  response 
amplitude  operators  can  be  used  In  such 
predictions  just  as  well  as  those 
obtained  from  regular  wave  tests. 

Professor  Tafcezawa  and  his  group  at 
Yokohama  National  University  (11,  12) 
are  continuing  to  develop  the  "tran¬ 
sient  wave  method"  discussed  In  another 
presentation  In  this  session  (13).  The 
transient  wave  method  has  the  advantage 
of  a  continuous  spectrum  In  both  the 
stationary  and  encountered  frequency 
domains.  Professor  Johnson  and  a  group 


from  both  the  Naval  Academy  and  MTS 
Systems  Corporations  have  developed 
the  "Periodic  Irregular  Encountered 
Wave  Technique"  (14)  which  Is  deter¬ 
ministic  (periodic)  to  the  moving  model, 
but  Is  stocastlc  (never  repeating)  at 
the  waveboard  and  In  the  tank.  This  Is 
accomplished  by  equally  spacing  the 
harmonic  components  In  the  encounter 
frequency  domain  so  that  the  Irregular 
wave  repeats  exactly  with  a  period 
equal  to  the  Inverse  of  the  frequency 
spacing.  The  on-line  computer  which 
drives  the  waveboard  calculates  the  non¬ 
equal  ly  spaced  component  frequencies 
which  are  then  summed  Into  a  never- 
repeating  drive  signal  for  the  wave 
generator.  A  different  drive  signal  Is 
required  for  each  carriage  speed  since 
the  frequency  shift  varies  with  the 
model  velocity: 
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IV. 


EXTREME  WAVE  AND  WAVEGROUP 
GENERATION  AND  ANALYSIS 


Another  area  covered  by  this  review 
Is  that  of  the  generation  and  analysis 
of  extreme  waves  and  "wave  groups". 

This  type  of  wave  action  Is  of  impor¬ 
tance  In  severe  wave  load  and  capsize 
tests  for  ship  models  and  for  the 
testing  of  breakwaters  and  other  coastal 
structures  In  hydraulic  laboratories. 

An  extensive  project,  "Ships  In  Rough 
Seas",  has  been  underway  at  the  Norwegian 
Hydrodynamic  Laboratories  since  1978 
(15-22).  The  project  Includes  detailed 
studies  of  breaking  waves  by  Kjeldsen 
and  Myrhaug  which  have  produced  new 
Insights  concerning  extreme  wave 
phenomenon . 

For  either  stationary  or  moving 
models  It  Is  frequently  desired  to  make 
a  wave  break  precisely  at  the  location 
of  the  structure  or  ship  model.  The 
methods  utilizing  sweep  frequency 
modulation  to  generate  transient  waves 
may  be  used  for  this  purpose  (18,  21), 
Including  those  for  repeating  transient 
waves  (13).  An  equivalent  method  for 
hydraulics  laboratories  has  been 
developed  by  Funke  and  Mansard  of  the 
NRC  Hydraulics  Laboratory  In  Canada  (24). 
Their  report  contains  a  copy  of  the 
SPLISH  software  with  sufficient 
remarks  to  reproduce  It  In  other 
laboratories.  Other  useful  wave 
generation  programs  are  also  available 
from  NRC  and  are  summarized  In  (25). 

A  novel  means  of  generating  extreme 


waves  by  using  a  traveling  flap  fitted 
to  a  towing  carriage  was  reported  by 
Hogben  (26). 

Extreme  waves  may  also  be  generated 
at  a  particular  location  through  the  use 
of  Inverse  Fourier  Transform  techniques 
with  predetermined  phase  control  rather 
than  randomizing  the  phases. *  In  the 
simplest  form  of  Inverse  FFT  generation, 
the  spectral  densities  of  each  component 
are  specified  and  the  group  velocity  of 
each  of  the  components  Is  computed  In 
order  to  predict  the  appropriate  Initial 
phase  relationships  to  produce  an 
additive  super-position  of  all  the 
frequency  components  at  a  particular 
location  In  the  wave  basin. 

Wave  grouping  Is  also  possible 
using  the  frequency  domain  developed 
for  the  Periodic  Irregular  Encountered 
Wave  Technique  (9).  In  this  technique, 
one  has  the  option  of:  a)  equally 
spaced  harmonic  frequencies  can  be  used 
to  produce  a  breaking  wave  at  a  parti¬ 
cular  location  every  repeat  period,  or 
b)  unequally  spaced  component  fre¬ 
quencies  can  be  used  to  produce  a  never 
repeating  wave  which  breaks  just  once 
at  the  desired  location.  If  the  equal 
frequency  spacings  are  chosen  in  the 
encountered  frequency  domain,  the 
program  can  produce  a  breaking  wave  on 
the  moving  model  once  each  repeat 
period.  In  the  latter  case,  the  extreme 
wave  location  moves  toward  the  wave- 
maker  at  the  speed  of  the  carriage  and 
does  not  occur  anywhere  else  In  the 
basin.  The  program  can  start  the  high 
frequencies  first  and  systematically 
add  in  the  lower  frequencies  to  avoid 
beach  reflections  as  long  as  possible. 

At  this  point,  a  distinction  should 
be  made  between  the  types  of  breaking 
waves.  The  largest  wave  forces  and  the 
greatest  capsize  dangers  for  smaller 
vessels  are  associated  with  plunging 
breakers  which  are  readily  produced  by 
shoaling  the  wave  Into  shallow  water 
(28)  or  by  narrowing  the  vertical  sides 
of  a  wave  flume  (29).  Wave-wave  Inter¬ 
actions  Involving  very  steep  Stokes 
waves  can  also  produce  plunging  breakers 
In  deep  water  as  reported  by  Kjeldsen, 
et.al.  (18-22).  Spilling  breakers  and 
deep  water  bores  (defined  In  reference 
18)  are  more  commonly  found  for 


♦Footnote:  Interestingly,  the  proper 
choice  of  random  phases  to  Insure  that 
breaking  waves  would  not  occur  for 
standard  ship  motion  tests  was  discussed 
by  Anderson  and  Johnson  (27).  To  deter¬ 
mine  a  response  amplitude  operator,  the 
spectrum  should  be  essentially  constant 
throughout  the  tank.  Breaking  waves 
distort  the  achieved  spectrum  down¬ 
stream  of  the  break  point. 


deterministically  generated  Irregular 
waves  In  most  deep  water  basins. 

However,  a  current  created  In  the 
opposite  direction  to  the  wave  direction 
will  steepen  the  wave  sufficiently  to 
cause  It  to  MplungeM  (18,21). 

Although  wave  spectra,  significant 
wave  heights  and  model  periods  are 
generally  reported  for  both  ocean  data 
and  for  model  testing  In  waves,  these 
parameters  are  not  adequate  to  properly 
describe  tests  Involving  extreme  waves 
or  "wave  groups".  The  averaging 
process  tends  to  mask  unusual  "events" 
and  several  radically  different  time 
series  can  all  be  represented  by  the 
same  spectrum  (30).  Spectral  shape 
parameters  are  also  Inadequate  for 
describing  events  which  are  related  to 
the  effects  of  phase  on  the  superposi¬ 
tion  of  wave  harmonics.  An  adequate 
characterization  of  "episodic"  waves 
and/or  "groupi ness"  appears  to  require 
a  time  domain  analysis  of  the  events. 
Three  approaches  to  such  an  analysis 
have  recently  been  put  forward.  The 
first  method  was  proposed  by  Kjeldsen 
and  Myrhaug  In  1978  (15,17,21).  It 
involves  the  use  of  the  "zero  down- 
cross"  method  of  wavehelght  analysis 
and  the  probability  distributions  of 
four  wave  shape  parameters.  They  are 
the  crest  front  steepness  parameter,  c  • 
the  crest  rear  steepness  parameters,  5  , 
the  vertical  asymmetry  factor,  X  ,  and 
the  horizontal  asymmetry  factor,  y  . 

They  analyzed  25,000  waves  and  report 
the  results  of  their  statistical 
analysis.  It  Is  not  yet  apparent  how 
to  reduce  these  characterizations  to 
simple  numbers  for  the  purpose  of 
specifying  a  desired  irregular  wave 
record. 

The  second  method  Is  the  "Smoothed 
Instantaneous  Have  Energy  History"  or 
SIWEH  as  proposed  by  Funke  and  Mansard 
(10)  In  1979.  From  the  SIWEH,  which  has 
a  time  axis  the  same  length  as  the  wave 
record,  Funke  and  Mansard  have  proposed 
a  Grouplness  Factor  which  equals  the 
standard  deviation  of  the  SIWEH  about 
Its  mean,  normalized  by  the  area  under 
the  wave  spectrum.  They  also  propose 
that  three  other  parameters  related  to 
wave  grouping  be  specified:  the  average 
group  length,  average  group  repetition 
period,  and  some  measure  of  the  group 
repetition  period.  This  Is  a  develop¬ 
mental  proposal,  and  operational 
definitions  of  all  the  parameters  have 
not  yet  been  worked  out. 

The  third  approach  to  a  time  domain 
analysis  of  the  extreme  wave  phenomena 
Is  the  "Half-Cycle  Counting  Technique" 
proposed  by  Buckley  at  the  19th 
American  Towing  Tank  Conference  (31). 
This  technique  was  adapted  from  a 
computerized  procedure  to  accumulate 
fatigue  load  data  In  structural  tests. 


The  wave  record  Is  banded  Into  uniform 
amplitude  Intervals  on  either  side  of 
the  mean  value  or  still  water  level. 

The  peaks  (both  maxima  and  minima,  thus 
the  term  half  cycle)  are  then  detected 
within  each  elevation  Interval  and 
summed  Into  the  appropriate  bln 
depending  on  whether  the  data  peak  Is  a 
"first"  peak  or  a  "second"  peak.  The 
counts  stored  In  the  bins  form  a  "data 
matrix"  which  Is  presented  In  a  visually 
satisfying  format  that  simultaneously 
logs  the  "half  cycle  counts",  the 
amplitude  histograms  of  the  occurrences, 
the  peak  count  occurrences,  and  any 
"episodic"  events  which  lie  outside  the 
expected  range  of  wave  elevations  for  a 
particular  sea  state.  This  procedure  is 
also  developmental  at  the  present  time 
and  has  the  choice  of  the  proper  dead 
zone  band  (threshold)  as  an  unresolved 
problem.  Neverless,  It  has  generated 
quite  a  bit  of  Interest  from  those 
exposed  to  Its  Interesting  presentation 
of  Irregular  wave  data. 

What  Is  currently  needed  Is  a  data 
analysis  package  which  Included  both 
the  Grouplness  Factor  and  Half-Cycle 
counting  technique  for  analyzing  wave 
groups  and  extreme  waves.  More  analysis 
of  actual  wave  records  Is  required  In 
order  to  define  appropriate  wave 
grouplness  parameters.  Representative 
wave  data  which  is  known  to  contain 
extreme  waves  and  wave  groups  should 
be  analyzed  by  various  tanks  to  get  a 
consensus  on  the  usefulness  and  adequacy 
of  these  approaches  before  any  standards 
are  recommended. 


V.  SHORT  CRESTED  IRREGULAR  WAVE 

GENERATION  AND  ANALYSIS 

A  state-of-the-art  summary  of 
multi-directional  wave  generation  and 
analysis  was  presented  at  the  15th  ITTC 
by  Hogben  (32).  Recent  progress  in 
methods  for  the  generation  of  short 
crested  seas  In  a  wave  basin  through  the 
use  of  multiple  wave  generators  1$ 
summarized  In  a  contribution  to  the 
16th  ITTC  Seakeeping  Committee  by 
Naeser  (33).  Three  methods  are  dis¬ 
cussed  for  producing  directional  spectra. 
The  "snake  principle",  pioneered  by 
NSMB,  has  been  expanded  to  multi¬ 
frequency  Inputs  at  the  University  of 
Edinburgh.  Their  wave  tank  Is  equipped 
with  89  wavemakers  which  may  be 
Independently  driven  through  the  use  of 
a  minicomputer  (34).  The  "Wallingford 
principle"  which  utilized  ten  Indepen¬ 
dently  controlled  wavemakers  arranged 
In  a  crescent  to  produce  a  desired 
directional  wave  spectrum  at  the  center 
of  the  basin  Is  described  by  Huntington 
and  Thompson  (35).  The  "diffraction 
principle"  also  uses  Independently 
driven,  snake  type  wavemakers, 
arranged  along  one  wall  In  combination 


with  reflecting  side  walls  rather  than 
beaches  (36).  Naeser  reports  that  this 
principle  will  be  tested  using  the  144 
wavemakers  arranged  along  the  63  m  wall 
of  the  Ocean  Laboratory  in  Trondheim 
In  1981  (33). 

Naeser  also  points  out  that  none 
of  the  techniques  is  able  to  produce 
a  mul ti -di rectional  spectrum  which  Is 
continuous  In  both  frequency  and  direc¬ 
tion.  The  snake  method  produces  a 
discontinuous  spectrum  containing 
discrete  frequencies  and  discrete 
di rections. 

The  Wallingford  method  produces  a 
nearly  continuous  frequency  spectrum 
containing  discrete  directions,  while 
the  diffraction  method  produces 
discrete  frequencies  from  "quasi 
continuous"  directions. 

The  analysis  of  mul tl-dl rectfona 1 
spectra  is  not  a  simple  task.  The  most 
recent  discussion  of  this  problem  Is 
that  of  Mathelsen  (37)  which  presents 
the  theory  and  a  computer  program  for 
the  calculation  of  directional  wave 
spectra  from  measurements  of  wave 
heights  at  the  corners  of  a  square.  The 
well  known  review  of  techniques  for 
directional  wave  spectra  by  Panicker 
(38)  1$  recommended  reading  for  those 
involved  In  such  studies. 
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DISCUSSION  TO  INVITED  PAPER  BY 
BRUCE  JOHNSON 


RESPONSE  TO  THE  DISCUSSIONS  CONCERNING  THE 
STATE  OF  THE  ART  REVIEW  OF  IRREGULAR  WAVE 
GENERATION  AND  ANALYSIS 


-  By  Stephen  J.  Rowe 

National  Maritime  Institute,  England  -  Bruce  Johnson 


Mr.  Chairman  I  have  two  points. 

Firstly,  I  am  concerned  about  the 
concept  of  the  generation  of  waves  using 
the  amplitude  modulation  of  the  component 
sine  waves,  and  the  confusion  this  intro¬ 
duces.  It  seems  to  me  that  all  this 
achieves  is  the  extension  of  the  time  series 
before  repetition,  (the  length  of  the  wave 
sequence  becomes  a  complex  function  of  the 
modulation  frequencies).  If  one  were  to 
record  such  a  time  series  and  perform  a 
Fourier  transformation  and  obtain  the 
coefficients,  one  could  in  principle  set  up 
a  mu  1 ti -component  sine  wave  synthesizer  or  a 
shift  register  synthesizer  in  amplitude  and 
phase  to  produce  exactly  the  same  time 
seri es .  The  amplitude  modulation  therefore 
is  effectively  simply  inserting  additional 
lines  into  the  spectrum,  and  one  will  of 
course  have  to  analyse  long  runs  in  the  tank 
to  make  use  of  the  finer  spectral  resolution 
that  has  been  thus  obtained. 

If  one  looks  at  the  time  history  of 
such  a  long  time  series  (however  produced) 
it  will  'appear  non-stationary'. 

Physically,  this  is  because  adjacent  lines, 
or  groups  of  lines  in  the  spectrum  are 
coming  to  anti-phase,  and  thus  cancelling 
out  the  energy  for  this  frequency  band  for 
a  period  of  time. 

Finally,  I  would  like  to  emphasize  that 
in  my  view  there  is  absolutely  no  point  in 
averaging  the  analysis  results  from  two  or 
more  repeating  sequences  of  waves.  All  this 
will  tell  you  is  the  repeatability  of  your 
wavemaker  and  instrumentation.  (You  may 
prefer  not  to  know  this!).  Worse  still,  if 
you  transform  as  one  block  two  or  more 
sequences,  your  analysis  resolution  will  be 
closer  than  the  'wet  lines'  in  the  tank  and 
you  will  over-resolve  the  spectrum.  It  is 
essential  therefore  that  one  should  generate 
long  sequences  of  waves  that  are  as  long 
as,  or  longer  than,  the  total  length  of  run 
record  that  you  wish  to  analyse. 


DISCUSSIONS  TO  PAPERS  BY  PROFS.  JOHNSON  AND 
TAKEZAWA 

-  By  Hun  Choi  Kim 

Ship  Research  Station,  KIMM,  Korea 

Irregular  wave  generation  assumes  an 
applicability  of  linear  superposition,. 

There  would  be  non-linearity  arising  from 
wave  decay  along  the  length  of  a  towing  tank, 
decay  being  dependent  not  only  on  each  wave 
component  but  on  the  whole  system  of  waves. 

Particularly  once  the  wave  or  the  wave 
system  breaks  there  would  occur,  it  would 
seem,  almost  certain  reduction  in  amplitude(s) 
as  well  as  phase  shifts  in  certain  com¬ 
ponents. 

I  wonder  if  such  has  been  noted,  and  if 
so,  how  serious  it  is  and  how  one  should 
go  about  in  compensating. 


Mr.  Chairman,  Mr.  Rowe  has  raised  two 
very  interesting  points  which  I  have  had  to 
research  after  returning  home  from  the  ITTC. 
With  regard  to  the  first  point,  it  is  not 
my  understanding  that  any  time  series  can  be 
reproduced  by  use  of  a  Tour i er  transform 
and  a  mul ti -component  sine  wave  synthesizer 
or  shift  register  synthesizer.  Linear  wave 
phenomena  represented  by  a  periodic  time 
series  (however  long  the  repetition  period 
so  long  as  it  is  finite)  can  be  reproduced 
by  the  sum  of  sinusoids  technique.  Transient 
phenomena  which  begin  and  end  with  zero 
amplitude  will  have  a  valid  Fourier  transform, 
but  the  frequency  resolution  will  be 
infinitesimal  and  these  time  series 
require  an  infinite  number  of  component 
waves  to  be  reproduced  exactly.  All  other 
time  series  cannot  be  reproduced  exactly  - 
only  an  approximation  to  the  original 
series  can  be  produced.  If  the  wave 
generation  system  is  a  linear  system  (true 
for  small  amplitudes)  this  approximation 
will  be  quite  good.  Since  the  discussion 
concerned  large  amplitude  waves  which  may 
be  quite  non-linear,  the  reproduction  of  the 
original  time  series  will  not  be  possible 
using  Fourier  transform  techniques. 

Amplitude  modulation  of  the  individual 
components  does  more  than  "effectively 
inserting  additional  lines  into  the  spectrum". 
In  recent  attempts  to  generate  extreme 
waves,  I  have  been  unable  to  prevent 
premature  breaking  when  using  the  sum  of 
sinusoids  technique  with  constant  component 
amplitudes  to  try  to  obtain  phase  convergence 
at  a  particular  location  in  the  tank. 

Using  the  sweep  frequency  method  (transient 
wave  method)  in  which  the  component 
amplitudes  vary  from  zero  to  whatever 
amplitude  is  generated  at  a  particular 
frequency  during  the  sweep,  very  large 
plunging  breakers  can  be  produced  in  deep 
water.  Neither  of  these  generation 
techniques  approximate  a  real  storm  sea, 
which  must  contain  time  varying  component 
amplitudes  to  realize  observed  extreme 
waves.  I  also  question  whether  the  extreme 
values  predicted  for  a  stationary  Gausian 
random  process  are  adequate  for  predicting 
the  occurance  of  extreme  waves  in  non- 
Stationary  storm  seas.  Incidently,  I  am 
not  aware  of  how  one  tells  whether  or  not 
a  time  series  "appears  non-stationary"  just 
by  looking  at  the  time  history. 

With  regard  to  the  second  point,  I  agree 
that  averaging  the  analysis  results  from  two 
or  more  repeating  sequences  of  waves  has  no 
point  if  the  system  is  linear.  If,  however, 
the  system  has  a  non-linear  response  which 
may  involve  transient  decay  or  amplification, 
the  analysis  of  more  than  one  repeating 
sequence  will  establish  whether  or  not  a 
steady  state  response  has  been  established. 

I  also  discovered  another  reference  (1) 
to  random  testinq  techniques  which  should 
have  been  included  in  the  original  paper. 
Although  the  reference  addresses  the  testinq 
of  automotive  structures,  its  conclusions 
are  also  applicable  to  the  testing  of 
marine  structures.  The  authors  of  this 


paper  state  that  pseudo  random  techniques 
involving  sufficiently  long  repetition 
periods  to  insure  adequate  spectral 
resolution  of  the  phenomena  being  observed 
are  the  most  efficient  for  linear  systems. 
However,  pure  random  testing  is  considered 
to  be  better  for  non-linear  systems  since 
"it  gives  the  best  linear  approximation  of 
a  non-linear  system."  Impact  testing 
(transient  method)  gave  good  results  for 
linear  systems,  but  is  less  effective  for 
non-linear  systems,  unless  only  a  "quick 
look"  estimate  is  desired. 

Or.  Schmiechen  has  raised  a  valid 
point  about  the  use  of  mathematically 
defined  time  series  analysis  terms.  Since 
there  are  semantic  difficulties  arising 
from  the  use  of  these  terms  to  describe  and 
categorize  wave  time  histories,  the 
Information  Committee  should  consider  the 
inclusion  of  operational  definitions  of 
time  series  analysis  terms  in  the  updated 
Dictionary  of  Ship  Hydrodynamics. 

Dr.  Kim  discussed  the  influence  of 
wave  breaking  on  the  wave  system  at  various 
locations  in  the  tank.  This  problem  was 
mentioned  in  the  footnote  on  the  fourth 
page  of  the  paper,  but  the  author  is  not 
aware  of  any  means  of  compensating  for  the 
wave  breaking  once  it  has  occured. 
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